A TWO- VARIABLE INTERLACE POLYNOMIAL 



RICHARD ARRATIA, BELA BOLLOBASt, AND GREGORY B. SORKIN 

Abstract. We introduce a new graph polynomial in two variables. This "interlace" poly- 
nomial can be computed in two very different ways. The first is an expansion analogous to 
the state space expansion of the Tutte polynomial; the significant differences are that our 
expansion is over vertex rather than edge subsets, and the rank and nullity employed are 
those of an adjacency matrix rather than an incidence matrix. 

The second computation is by a three-term reduction formula involving a graph pivot; 
the pivot arose previously in the study of interlacement and Euler circuits in four-regular 
graphs. 

We consider a few properties and specializations of the two- variable interlace polynomial. 
One specialization, the "vertex-nullity interlace polynomial" , is the single- variable interlace 
graph polynomial we studied previously, closely related to the Tutte-Martin polynomial 
on isotropic systems previously considered by Bouchet. Another, the "vertex-rank interlace 
polynomial" , is equally interesting. Yet another specialization of the two- variable polynomial 
is the independent-set polynomial. 



1. The INTERLACE POLYNOMIAL 

In |A HSnOal EHSj, we introduced a single-variable "interlace" graph polynomial. It 
emerged that the interlace polynomial could be regarded as a special case of the Tutte- 
Martin polynomial of an isotropic system, as discussed briefly in Section |3] here and more 
fully in jABSj . We defined the polynomial by a recurrence relation, and Balister, Bollobas, 
Cutler and Pebody |BBCP02'] used a property of the nullities of the adjacency matrices 
of the graphs in question to resolve a conjecture posed in [ABSOOa'. The hnear-algebraic 
approach of fB BCPOS] . extended to embrace the matrix ranks as well as nullities, led us to 
the two- variable polynomial introduced here. The two-variable interlace polynomial appears 
to be something entirely new; in particular, we are aware of no two-variable versions of the 
Tutte-Martin polynomial or other closely related polynomials. The interlace polynomial is 
an entirely different object from the Tutte polynomial, but the two have an extremely inter- 
esting structural similarity which immediately suggests a family of additional polynomials 
for further exploration. We will now define the two-variable interlace polynomial. 

Given a graph G with vertex set V{G), for any subset 5* C V{G), let G[S] be the subgraph 
of G induced by S. We allow graphs with loops on their vertices, and we also allow the null 
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graph with no vertices, but we do not allow multiple loops or multiple edges. Write Q for 
the set of graphs including the null graph. 

For a matrix A over F2, let n{A) be the nullity of A and r{A) its rank. Abusing notation 
slightly, for a graph G, n{G) and r{G) will denote the nullity and rank of its adjacency 
matrix, so n{G) + r{G) = |V^(G)|. (For example, if S is the empty set then G[S] is the null 
graph of rank and nullity 0.) We remark that for loopless graphs G, r{G) is always even, as 
the rank of a zero-diagonal symmetric matrix. 

We define the two-variable interlace polynomial q{G; x, y) of a graph G of order n as a 
sum of 2" terms: 

q{G;x,y)= (x - l^^^'Hy - 1^""^'^ (1) 

SCV(G) 

the sum taken over all subsets including 5 = and S = V{G). For convenience we define 
the monomial 

m{H) = 

so that 

g(G) = 5^m(G[5]). (2) 

s 

This "state space" expansion of the two- variable interlace polynomial may be seen as an 
analogue of the Tutte polynomial given as 

T(G; x,y) = J2ix- - l)"^^). (3) 

FCE 

The significant differences are that our sum is over induced subgraphs of G (as given by vertex 
subsets) rather than arbitrary subgraphs on the full vertex set (as given by edge subsets), and 
our rank and nullity are those of the subgraph's adjacency matrix rather than its incidence 
matrix. (In the context of the Tutte polynomial, a graph's rank is normally defined as the 
number of vertices minus the number of components, and it is easy to check that this is the 
F2-rank of the incidence matrix.) That our subgraph rank appears positively rather than 
subtracted from the rank of the whole is not significant, as it can be adjusted by a change 
of variables. That is, if one prefers the polynomial Y1scv{G)(-'^ ~ lY'''^^~'^^'^^^\y — 1)"'('-^I'^1), 
it is just (x-l)'^(«)g(^,|/). 

A surprising basic property of this polynomial is that for loopless graphs it satisfies a 
three-term reduction formula, as per Theorem |2l and for looped graphs, a pair of reductions, 
per Theorem El we prove these results in Section |21 In Section [HI we show a pair of properties 
of the two- variable interlace polynomial. In Section |3] we describe the polynomial's special- 
izations to the vertex-rank polynomial and the vertex-nullity polynomial, and in Section 
its specialization to the independent-set polynomial. We calculate the polynomial on some 
basic graphs in Section IHl We conclude in Sections [7| and |H1 with thoughts on generaliza- 
tions of the polynomial, some of which seem likely to prove interesting, and with some open 
problems. 
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We now proceed to the interlace polynomial's reduction formula. 

2. Reduction formula 

We begin by showing a reduction involving a "pivoting" operation on an edge ab of G for 
which both a and b are loopless. It is natural to start with this case because it provides 
a recursive definition of q{G) for the natural class of loopless graphs G. In the following 
subsection, we show how to reduce on a looped vertex a of G, completing a recursive definition 
of q{G) for arbitrary graphs. 

2.1. The pivot and reduction. As in |ABSnnbl IXHRl IBou99j . and related to Kotzig's 
transformations on Euler tours |Kot69j , for a graph G and an ordered pair ab = (a, b) of 
distinct vertices of G, we define the pivot operation on ab mapping G into G"^ as follows. 
We say that two vertices x,y of G are distinguished by {a, b} if x,y ^ {a, b} and x, y have 
distinct non-empty neighborhoods in {a,b}. Let G"* be the graph with vertex set V{G) in 
which xy is an edge if either xy ^ E{G) and x and y are distinguished by {a,b}, or else 
xy e E{G) and x and y are not distinguished by {a, b}. 

Let us spell out this definition in detail. Partition the vertices other than a and b into 
four classes: 

CI: vertices adjacent to both a and b; 
C2: vertices adjacent to a alone; 
C3: vertices adjacent to b alone; and 
C4: vertices adjacent to neither a nor b. 

Definition 1 (Pivot). A graph G is pivoted on vertices a, b to obtain G""^ as follows. For any 
vertex pair xy where x is in one of the classes (C1)-(C3) andy is in a different class (Cl)- 
(C3), the pair xy is "toggled": if it is an edge of G it is not an edge of G""^ , and if it is not 
an edge of G then it is an edge of G^^ . All other pairs of vertices are adjacent in G""^ iff they 
are adjacent in G. 

Trivially, pivoting and restriction to a subgraph satisfy a commutative law: for any S C 
V{G) with a, 6 G ^, 

G''\S] = {G[S]yK (4) 

Although pivoting is defined for any vertex pair ab, we shall only exploit it in cases where 
ab is an edge and a and b are both loopless. 

We shall write the adjacency matrices of G and G"'^ with rows and columns put into six 
groups according to their relations to a and b. The first group consists of a alone, and the 
second of b alone; groups three to six are the four classes above. Write 1 for an all-1 row or 
column vector of whatever dimension and likewise for an all-0 vector. Then, for a graph 
G with loopless vertices a and b and having an edge ab, the adjacency matrix of G is of the 
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form 
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(5) 



where in all cases Mjj is the transpose of Mjj. Then for such a graph the adjacency matrix 
of G''^ is 



A{G 
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(6) 



where M'^ denotes the complement of the matrix M. 



-h). 



Lemma 2. For any graph G with an edge ah, with a and h both loopless, r{G—a) = r{G' 
andr{G) = r{G°'^ — a — h) + 2; equivalently, n{G — a) = n{G'^^ — a) andn{G) = n{G°'^ — i 

Proof. With the adjacency matrix A = A{G) represented as in (jSj), add row 2 to each row 
in the 3rd and 4th groups. Repeat the same operations on columns to give a matrix 
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Note that, outside of the neighborhoods of a and h (the first two rows and columns), A' is the 
adjacency matrix of G"^; moreover, since these linear operations are invertible, r[A) = r{A'). 

To prove the first assertion, discard the first row and column of A to yield A\a and 
similarly that of A' to obtain A'\a. Since they did not use row or column 1 (vertex a), the 
same linear transformations as before map A\a to A'\a, showing that r{A\a) = r{A'\a). 



But A\a = A{G-a) and A'\a = A(G"^^ - a), therefore r{G 



r{G' 



ab 
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To prove the second assertion, further transform A' by adding row 1 to each row in the 
3rd and 5th groups, and repeating for columns, to obtain 



A" = 



( ° 
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Because these hnear transformations are all invertible, they preserve rank and nullity: r{A) = 
r{A"). Referring to ©, note that A"\{a,b} is the adjacency matrix of G""^ — a — b. The 
first and second rows in A" are linearly independent of one another and of all other rows, so 
deleting them reduces the rank by 2. After deletion of these rows the first two columns are 
all-zero, so deleting them does not change the rank; it follows that r{A") = r{A"\{a, b}) + 2. 
That is, r{G) = r{A) = r{A") = r^C^ - a - 6) + 2. □ 

Note that if either a or 6 has a loop, the top-left submatrix of A{G) differs from that in (0), 
resulting in a different "border" in the matrix A", so that the border's deletion changes the 
nullity unpredictably. That is, if there is a loop at either a or b, n{G) — n{G°'^ — a — b) may 
be rather than 2. Thus, the looped case is dealt with in Section IT^ 

Theorem 3. For any edge ab of a graph G, where a and b are both loopless, 

q{G) = q{G - a) + qiG"^ - b) + ((x - 1)^ - l)q{G^^ -a-b). 

Proof. For 5* ranging over subsets of V{G)\{a,b}, by (0), 

q[G) = ^{m(G'[5]) + m{G[S U a]) + m{G[S U 6]) + m{G[S U {a, b}])} (7) 

s 

while 

q{G -a) + q{G^^ -b) + ((x - 1)^ - l)g(G'"* -a-b) 

= ^{m((G - a)[S]) + m((G - a)[S U b])} 
s 

+ ^{m(G"^^ - b)[S]) + m((G'"^ - b)[S U a])} 
s 

5 

To show that ((Tj) and (jHI) are equal, we will show equality of their terms for each S. Two terms 
of ((Tj) directly match their counterparts in (jH)): m{G[S]) = m{{G — a)[S]) and m{G[S Ub]) = 
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■m{{G — a)[S U b]) . A third equality follows from the first part of Lemma |21 

m{G[S U a]) = m{G[S U {a, b}] - b) 

= m{{G[SVJ{a,b}]T^ -b) 
= m^G^^S a]) (by©) 
= m((G'^^ -6)[SUa]). 

The final equality, between a single term from (jZj) and two terms from (jHl), follows from the 
second part of Lemma |21 

m{G[S U {a, b}]) = {x- lfm{{G[S U {a, b}]Y^ - a - b) 

= {x- l)='m(G'"^[5]) 

= m(G"^'[5]) + {{x - 1)2 - l)m(G°'[^]) 

= miiG"'' - b)[S]) + {{x - 1)2 - l)m((G'"' -a- b)[S]). 

□ 

2.2. Local complementation. In the case graphs with loops, there may not always be an 
edge ab, with both a and b loopless, to which Theorem El (derived from Lemma 12)) may be 
applied. In this case, though, there must be a looped vertex, and we may apply a different 
reduction instead. 

Bouchet |Bou99j defines the "local complement" G"^ of graph G on vertex a by comple- 
menting (toggling the presence or absence of all edges, including loops) the subgraph of G 
induced by the neighborhood of a, while keeping the graph otherwise unchanged. To be 
precise, let T{a) be the set of neighbors of a; in particular a G r(a) iff there is a loop on a. 

Definition 4 (Local complementation) . A graph G is locally complemented on a vertex a to 
yield G"" , where is equal to G except that G'"[r(a)] = G'[r(a)]. 

In notation like that of (0) and © but only distinguishing a vertex a having a loop, its 
neighbors, and its non-neighbors, we may write 
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1 





1— ' 


Mn 


M12 





M21 


M22 



A{G) = I 1 Mn M12 I (9) 

and 
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h- ' 
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Mfi 


M12 





M21 


M22 



A{G^) = 1 Mfi M12 . (10) 
\ M21 M22 / 

Incidentally, it is observed in |Bou99j that a pivot is equal to a composition of local 
complementations, G""^ = {{G°'Y)°', followed by a swap of the labels a and b. (This is correct 
as stated for our version of pivoting, which differs from Bouchet's by a label swap.) 
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Lemma 5. For any graph G with a looped vertex a, r{G) = r^G"" — a) + 1, and equivalently 
n{G) =n(G"-a). 

Proof. With the adjacency matrix A = A{G) represented as in (jHl), let A' be obtained by 
adding row 1 of A{G) to each row in the second class, and repeating for columns; thus 

/ 1 

A' = Mfi Mi2 
\ M21 M22 

Referring to note that A{G°' — a) = A{G°')\a = A'\a. The linear transformations are 
invertible, so r{A') = r{A). The first row of A' is independent of the others, so deleting it 
decreases the rank by 1, and what remains of the first column is all-zero, so deleting it does 
not change the rank; thus r{A'\a) = r{A') — 1. We conclude that r(G"' — a) = r{A'\a) = 
r{A')-l = r{A)-l=r{G)-l. □ 

Theorem 6. For a graph G, for any edge ah where neither a nor h has a loop, 

q{G) = q{G -a) + g(G"^^ - b) + ((x - 1)^ - l)g(G"^'' -a-b), (11) 

and for any looped vertex a, 

q{G) = q{G-a) + {x-l)q{G''-a), (12) 

Proof. Equation (fTTj) is precisely Theorem El repeated here to give one comprehensive theo- 
rem. 

To prove (fT^ . as in the proof of Theorem El we show that for each S C V{G)\a, we have 
equality between the corresponding summands of 

q{G) = ^{m(G[^]) + m{G[S U a])} 
s 

and those of 

q{G -a) + {x- l)m(G'" - a) = ^m((G - a)[S]) + {x - 1) ^m((G" - a)[S]). 

s s 
The first terms, m{G[S]) and m((G — a)[5']), are identical. By Lemma m{G[S U a]) = 
{x — 1) m{{G[S U a]y — a) = (x — 1) m{{G°' — a)[S]), which completes the proof. □ 

The reduction formulas give an alternative characterization of the two-variable interlace 
polynomial. We write En for the empty graph on n vertices. 

Corollary 7. The two-variable interlace polynomial defined by (0) is the unique map q : Q ^ 
Z[x, that satisfies the reduction formulas illil^) and the boundary conditions q{En) = y"", 
n = 0,l,.... 

Proof. By Theorem El the two- variable interlace polynomial q{G) defined by (Q) satisfies the 
reduction formulas (jllll2p and from ([Q) it is immediate that it also satisfies the boundary 
conditions. Uniqueness follows from pill2j) by induction on the order of the graph. □ 
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3. Two PROPERTIES 

As will be shown in the next section, the single-variable interlace polynomial defined 
m [XBSOObJ. which there was denoted as x), is a special case of the present two- variable 
interlace polynomial, and we will denote it here as qN{G] y). (The reason for this notation, 
and the definition of ^at, will be given in the next section.) 

Since |ABS00b] showed that this single-variable interlace polynomial satisfies the identity 
gAr(G; y) = qj\f{G"'^; y) (at least for loopless graphs, all that jABSOObj considered), something 
of the same sort might be expected for the two- variable polynomial. In fact it is not generally 
true that q{G] x, y) is equal to q{G"'^; x,y) — a counterexample is the path of length 3, pivoted 
on the middle edge — but instead we have the following proposition. 

Proposition 8. For any graph G with edge ah, 

q{G -a) -q{G-a-b) = q^G^"^ - a) - qi.G"^ -a-b). 

Proof. By the interlace polynomial's definition, with sums taken over subsets S C V{G)\{a, b}, 

q{G ~a) -q{G-a-b) = ^(m(G[5]) + m{G[S U b])) - ^m(G[5]) 

s s 

= J2MG[SUb]) 

s 

= ^m(G'"*[^U6]) 

s 

by Lemma |21 From this point symmetry completes the proof. □ 

As with the earlier single- variable polynomial qjs[{G;y), the present two- variable polyno- 
mial obeys a simple product rule. For graphs Gi = (Vi, Ei) and G2 = (V2, E2) with disjoint 
vertex sets, Vi fl V2 = 0, let Gi U G2 denote their disjoint union, (Vi U V2, -Ei U £'2)- 

Proposition 9. If Gi and G2 are graphs on disjoint vertex sets, and Gi U G2 their disjoint 
union, then q{Gi U G2) = q{Gi)q{G2) ■ 

Proof. As in |ABSnnbj . the relat ion can be proved through the reduction, taking advantage 
of the fact that pivots in one component do not affect the other. Alternatively, it is easy 
to see that the sum ((T)) is the product of the corresponding sums over all Si C V{Gi) and 
^2 C V{G2). □ 

4. Specializations of the interlace polynomial 

Specializing g(G; x, y) by setting x = 2 (or, for loopless graphs, where the rank is always 
even, x = 0) causes the "rank" term to disappear from (Q), and so we call this polynomial 
the "vertex-nullity interlace polynomial", 

q^{G-y) = q{G-2,y)= J] (y - 1)"(^H). (13) 

SQV{G) 
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As alluded to in the previous section, this is the single- variable interlace polynomial stud- 
ied in lAi^SOOl) ■ To see this, note that substituting x = 2 into (fTTj) gives the reduction 
qN{G]y) = qj\f{G — a;y) + qj^iG""^ — h;y), precisely the reduction that defined the single- 
variable polynomial in |ABS00b] (which was defined only for loopless graphs). The boundary 
conditions qj\f{En;y) = y"' also match, concluding the equivalence. 

Corollary 10. The single-variable interlace polynomial of |ABS00b] . defined there by a 
reduction formula, has the explicit expansion ()13|) . 

This single- variable interlace polynomial qN{G) is related to the Martin polynomial and 
circuit partition polynomials jMar781 IBou871 IBou88[ IBouQll IEM98L lEMnC). .Bol02l \KB^. 
Bouchet IBouOOj recognized that the single- variable interlace graph polynomial of |ABS00b] 
was a specialization of the Tutte-Martin polynomial of an isotropic system (introduced by 
Bouchet in |Bou87j and generalizing the Martin polynomial), and this connection was clari- 
fied and made explicit by Aigner and van der Hoist in |Avd Hn4] . 'Av dHn4] . written after the 
present work's submission, also proved a conjecture from ABSOOaJ that gAr(G; —1) is always 
of the form ±2^* (previously proved in [BB CP02] ). independently derived the expansion (fTSj) 
of the nullity polynomial, and introduced a related polynomial. However, in contrast to 
the Tutte-like two-variable graph polynomial of the present paper, we are not aware of any 
two-variable generalization of the Tutte-Martin polynomial. 

Analogously to the nullity polynomial qN{G] y), there is also a single- variable "vertex-rank 
interlace polynomial" 

qn{G-x)=q{G-x,2)= J] (x - 1)^(«[^1), 

SCV{G) 

and it appears to be equally interesting. First, the vertex-rank polynomial distinguishes 
graphs of small order better than the vertex-nullity polynomial. For example, the rank 
polynomial distinguishes all 11 simple graphs of order 4, where the nullity polynomial takes 
on only 8 distinct values. (And the rank polynomial distinguishes all 90 looped graphs of 
order 4, where the nullity polynomial takes on only 17 distinct values). At order 5 there 
are 34 non-isomorphic simple graphs: the rank polynomial takes 33 values, and the nullity 
polynomial only 17. (There are 544 non-isomorphic looped graphs, the rank polynomial takes 
541 values, and the nullity polynomial only 41.) Similarly for trees: the nullity polynomial 
fails to distinguish one pair of trees of order 8 and two pairs of order 9; the rank polynomial 
distinguishes all trees of orders 8 and 9. 

In HHE] it was shown that certain basic graph parameters could be read out from the 
vertex-nullity polynomial, namely the order, the component number, the edge-independence 
number, and an upper bound on the (vertex) independence number. The vertex-rank poly- 
nomial, too, gives the order. 

Remark 11. For any graph G of order n, qnlG; 2) = 2". 

Proof. The formula in (jT]) reduces to a sum, over all 2" subgraphs of G, of 1 raised to a 
power. □ 
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As per the following proposition, the maximum degree of either variable in the two- variable 
interlace polynomial is unchanged by "removing" the other variable (substituting 2). 

Proposition 12. For any graph G, deg^{q{G; x, y)) = deg(g/j(G; x)), and degy{q{G; x,y)) = 
deg(gAr(G; y)), where deg^, (respectively degy) denotes the maximum degree of x (resp. y) in 
the polynomial. 

Proof. We will prove the statement for the vertex-rank polynomial; that for the vertex-nullity 
polynomial is proved identically. Since qniG] x) = q{G; x, 2), deg{qR{G; x)) < deg^(g(G; x, y)) 
Consider any S C V{G) contributing to a term of the maximum x-degree, degree k. But 
qR{G;x) = J2scv{G)(^ ~ 1)''(<^['5'])^ so each such term here also has x-degree k. For each 
such S the coefficient of a;^ is 1: there is no cancellation, and so the x-degree is k in g^j as it 
was in q. □ 

In |ABSj we showed that deg(giv(G)) > ind(G), that is, deg(gAr(G)) is an upper bound 
on the independence number. (While it was proved there only for loopless graphs, the same 
result for looped graphs follows from (fT^ . in the next section.) We showed graphs for which 
the absolute gap deg(gAr(G')) — ind(G') was arbitrarily large, but it remained open whether 
the ratio deg(gAr(G'))/ ind(G) could be made arbitrarily large. Indeed it can. 

Let ifrf be the d-dimensional Hamming cube (so \iid\ = S'^) and let be its complement. 
Then we have the following. 

Remark 13. The Hamming cube and its complement have (respectively) clique and inde- 
pendence numbers cl{Hd) = ind{Hd) = 2. For d > 1 odd, their nullities are n{Hii) = and 
n(Hd) = 2'^/2, and for d even, n{Hd) = 2^/2 and n(Hd) = 0. 

For Hd with d > 1 odd, the set 5* = V{Hd) is a witness that deg{qN{Hd)) > 2'^/2, and 
therefore deg(giv(^d))/ md(Fd) > 2'^/4. 

Proof. That md{Hd) = cl{Hd) = 2 is immediate from the structure of the Hamming cube. 
We next derive the nullity of Hd, from which that of Hd follows quickly. Let the adjacency 
matrix of Hd be Ad and that of Hd be Ad- 

We first show that Ad is self-inverse, for odd d. The dot product of the ith and jth rows 
of Ad is the parity of the number of vertices at Hamming distance 1 to both i and j in Hd- 
If i = j, then there are d such vertices, for parity 1. If the distance between i and j is more 
than 2, then there are no such vertices: parity 0. And if i and j are at distance exactly 2, 
then the vertices in question must agree with both i and j where those two agree, agree with 
i in one of the two coordinates where it differs from j, and agree with j in the other; there 
are 2 such vertices, for parity 0. When d is even everything is the same except in the case 
i = j, where the d common neighbors mean parity 0, and thus = 0. 

The situation is similar for Ad, with the parities reversed: Modulo 2, Ad = Ad + 1 + I , 

2 2 

where 1 is the all-1 matrix. Then for d odd A^ = I, and for d even A^ = 0. 

For d odd, invertibility of Ad means it is of full rank. The self-invertibility of Ad for d 
odd can also be used to show that Ad is of half the full rank for d even: "Gluing" together 
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two copies of Hd-i to make an H^., = ( ^ } ) , self-invertibility of means 



the second set of rows is simply A^^x times the first set. Since the second set is a hnear 
combination of the first set, the rank of the whole matrix is at most 2"^"^; the presence of a 
block / means this rank is achieved. _ 

Similarly, for we have = I -r I ) where J = 1 + /. Here we find that for 

V Ad-\ J 

d odd, the second set of rows is just A^-iI = I A^-i times the first set, and conclude that 
r(Aa) = 2'^-\ □ 



5. Counting independent sets 



There are other interesting specializations of the 2- variable interlace polynomial. Evaluat- 
ing at y = 1 means that (y — = except when n(G[H]) = giving (y — l)"^^!^]) = 1. 
In particular, then, q{G; 2, 1) = Ehcg(2 - ^y^^Kl - 1)"^-^^ = Eh: n{H)=o l*"^^^ counts full- 
rank induced subgraphs of G. 

Similarly, q(G; 1, 2) = J2h- r{H)=o l"^^-* counts the independent sets of G (including the 
empty set), a problem that has received widespread attention. In particular, it is known 
that counting independent sets (computing the independence number) is T^P-complete even 
for low-degree graphs |DG00j . so it follows that it is #P-hard to compute the two- variable 
interlace polynomial (in particular at the point {x,y) = (1,2)) and the "rank" interlace 
polynomial (at a; = 1). In fact, it is hard to compute the independence number even 
approximately |DFJ98j . and so the interlace polynomial must be hard even to approximate 
at these points. 

Given the similarity to the Tutte polynomial, which is hard to compute almost everywhere 
( |J VW90] . see also |Wel93j for a survey), and given the variety of structures counted by the 
interlace polynomial (see |ABSj ) . with counting typically being #P-hard, it is anything 
but surprising that the interlace polynomial is computationally hard. However, it was a 
question left unresolved in |ABSj . and in fact we still do not have a proof that computing 
the "nullity" polynomial is #P-hard. Moreover, in analogy with the Tutte polynomial, it 
would be of interest to show that the interlace polynomial is hard to compute at almost all 
points (x, y). 

A particularly interesting evaluation is 

g(G; 1, 1 + A) = ^ O'-^^^A"^^) = ^ AI^I=/(G;A), (14) 

HCG H: r{H)=0 

the independent-set polynomial (the sum over all /c > of A'^ times the number of inde- 
pendent sets of cardinality k). It is well known that further quantities of interest can be 
computed from the independent-set polynomial and its derivatives. For example from (fT^ 
it is clear that ^q{G; 1, 1 + A) = E//: r(H)=o |-^^|A''^'"^ For A = 1 this is just E//: r{H)=Q l-^l' 
so that ^(1,2) is the sum of the sizes of all independent sets. 
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6. Polynomials of some basic graphs 

We compute the interlace polynomial of some basic graphs, notably complete graphs Kn-, 
complete bipartite graphs Km,n, and paths P„ of length n. 

Proposition 14. For all n and m we have 
q{En) = 

Q{Kn) = i (x" + (2 - xr) + \ [y—^^ (x" - (2 - xT) 

q{Km,n) = ((y™ - - 1)) + y'" + - 1 



1/ 3y + 2x(x-2) \ l + ^l+Aiy + xix-2)) 

2 [y ^l + 4{y + x{x-2))) [ 2 

^ 1 _ Sy + 2xix - 2) A / 1 - v/l + % + x(x-2))' 

Proof. That q{En) = is immediate from (^Q) and also figured into the boundary condition 
in Corollary [7| 
For Kn we have 

k<n ^ ^ 

which, letting odd(A;) = 1 if /c is odd and otherwise 

fc<n ^ ^ 



dd(fc) 



fc even ^ ' fc odd ' 



The even and odd sums are computable from the sum and difference of {z + 1)" = ^ 
and {—z + 1)" = substituting z = x — 1 and simplifying gives the expression 

shown. 

We also derive q{Km,n) directly from (P). Km,n has (7) ( subgraphs /^jj. Each such 
subgraph's adjacency matrix has the form 



1 

1 
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whose rank is 

Jo if i = or j = 
~ [2 ifi>Oandj>0. 

Then 

^Ep-m.-ir-(T)(;).E(;).E(T)-. 

the four terms coming respectively from the cases where i > and j > 0; i = and J > 0; 
j = and i > 0; and i = j = (the null subgraph). Expanding, 



(&-"(■;)) (s 

n ^ \ m 

+Efe-iK';)+E(i'-i: 

7=1 ^-^^ i=l 



+ 1. 



The claim for K^ n follows immediately. 

For G = Pn with n > 2, we use the reduction (jllj) with edge ab, where 6 is a leaf. 
Since G — a is the disjoint union of Pn-2 and Ei, q{G — a) = yq{Pn-2)- Since G"'^ = G, 
G"^ -h = Pn-i and G"'' - a - b = Pn-2- The net result is 

g(P„) = {y + x^- 2x)q{Pn-2) + q{Pn-i)- 

Solving this recursion, with the boundary conditions q{Po) = q{Ei) = y and q{Pi) = ^'(-^2) = 
x'^ — 2x + 2y, yields our formula for □ 

7. Further polynomials 

We observed in Section^ that the interlace polynomial's expansion is similar to that of the 
Tutte polynomial, with two significant differences: the sum is over vertex rather than edge 
subsets, and the rank and nullity are those of an adjacency matrix rather than an incidence 
matrix. 

This suggests a whole range of polynomials given by similar expansions, with the sums 
taken variously over vertex or edge subsets, with rank and nullity being those of an adja- 
cency or an incidence matrix, and with the field used perhaps being other than F2. Of course 
there is no obstacle to taking a "master polynomial" summing over both vertex and edge 
subsets, incorporating terms for both the adjacency-matrix and the incidence-matrix rank 
(using four variables instead of two), and even introducing further variables to incorporate 
ranks computed over other fields. It would be interesting to determine which of these poly- 
nomials satisfy reductions akin to that of Theorem El and which ones have significance in 
combinatorics or other fields. 
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8. Open problems 

As was the case with the single-variable interlace polynomial, the two-variable interlace 
polynomial is quite new, and there are more questions than answers. Here we simply list a 
few of the obvious ones. 

Is q{G) reconstructible, i.e., given q{G — a) for each vertex a, can we reconstruct q{G)l 

What is the expectation of q{G) for a random graph Gl 

We conjectured in |ABSj that the vertex-nullity polynomial's coefficient sequence might 
be unimodal. For the two-variable polynomial, representing the coefficients of —x, y) 
as an array whose entry is the coefficient of x^y^ , the array's rows and columns are 
also unimodal, for all (loopless) graphs through order 7. (With the substitution of —x for 
X, it is clear from Corollary [7| — but not from (Q) — that the coefficients are all non- 
negative.) Unfortunately, there are six graphs of order 8 for which this is not the case. Since 
however the vertex-rank polynomial and the vertex-nullity polynomial for these six graphs 
do have unimodal coefficient sequences, it remains possible that both of these single- variable 
polynomials' coefficient sequences are always unimodal. (As pointed out in |ABSj for the 



nullity polynomial, though, there are reasons to be doubtful, including the falsification in 
1993 of a similar conjecture for the Tutte polynomial |Sch93j .) 

But the most promising line of inquiry is the proposal in the previous section, to generalize 
from the interlace polynomial (Q) and the Tutte polynomial to generate other interesting 
polynomials. 
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